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Some recent measurements of the chemical composition of the cosmic radiation indicate that at the energy of 
3 x 10 18 eV, around the ankle, light cosmic ions dominate the spectrum as it occurs in the preknee energy region. 
Taking advantage of a recent theory of cosmic radiation which provides a quantitative explanation of the knee, the 
second knee and the ankle, the chemical composition of cosmic radiation is explicitly calculated giving individual 
ion spectra and ion fractions from 10 12 eV to 5 x 10 19 eV. The calculation assumes two components of the cosmic 
radiation feeding the ion flux at Earth: one originated in the disc volume and another one, called extradisc 
component, which from the disc boundaries traverses the Galaxy reaching the solar system. Data above 10 17 
collected during half century of experimentation by Auger, HiRes, Agasa, Akeno, Fly' s Eye, Yakutsk, Haverah 
Park and Volcano Ranch experiments are reviewed, examined and compared with the theoretical <ln(A)>. 
The comparison between computed and measured <ln(A)> exhibits a good global accord up to 2 x 10 19 eV except 
with the HiRes experiment and an excellent agreement in the range 10 15 -10 17 eV with Kascade, Eas-top, Tunka 
and other experiments. The accord requires a flux of the extradisc component of 1.8 x 10 14 particles/ (m 2 sr s 
eV 1 ' 5 ) at 10 19 eV, twice that generated by disc sources. 



1. Introduction 

The spectrum of the cosmic radiation be- 
tween the knee and the ankle has been recently 
calculated according to the Theory of Constant 
Spectral Indices [1-4]. The galactic mechanisms 
and the parameters generating the knee are the 
same at the origin of the ankle: this circumstance 
constitutes a major characteristic of the solution 
of the knee and ankle problem. Another notable 
aspect of this solution is the presence of a struc- 
ture at (5 — 7) x 10 17 eV which is called the sec- 
ond knee, just a direct consequence of the theory 
[5] . Figure 1 shows an example of excellent quan- 
titative accord between computed and measured 
cosmic-ray spectrum [4] . 

The Theory of Constant Spectral Indices de- 
termines the energy spectra of individual ions, or 
group of ions, in the interval 10 11 -5 x 10 eV 
with only one normalization point for the par- 



1 This paper is an english translation of an INFN Re- 
port to be published in italian: Composizione Chimica 
della Radiazione Cosmica intorno alia Caviglia e Indici 
Spettrali 



ticle flux. As a consequence, once the spectral 
indices and the abundances at a given energy are 
assigned, the chemical composition of the cosmic 
radiation is known at any energy. Fig. 2 shows 
that individual ion spectra in the knee energy re- 
gion are also in good agreement with the theory. 
The simplicity of the theory and its anchorage to 
the empirical parameters governing the cosmic- 
ion motion in the Galaxy are not minor aspects 
besides the excellent quantitative agreement with 
ion spectra exhibited in figures 1 and 2. 

The chemical composition expressed by 
<ln(A)> has been also calculated with all cosmic- 
ray sources placed in the disc [BJ, thus ignoring 
any extragalactic component. Figure 3 reports 
the results of the calculation. The profiles of 
<ln(A)> measured by the Kaskade and Tunka 
experiments along with that extracted with the 
measurements of X max of the Auger Collabora- 
tion are shown in figure 3. There is an evident 
disagreement above 10 17 eV between theory and 
data. All other experiments above 10 17 eV have 
been purposely omitted from figure 3 to render 
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Figure 1. The differential cosmic-ray spectrum 
(black thick line) resulting from the theory nor- 
malized to 3.79 x 10 17 particles/(m 2 sr s eV 15 ) 
at 10 12 eV (green square). The data are from the 
Kascade [7] and Haverah Park [8] experiments. 



the disagreement between theory and data more 
evident. According to the above quoted mea- 
surements the theoretical chemical composition, 
as generated by the galactic cosmic-ray sources 
above 4xl0 17 eV, is incompatible with the data, 
primarily because it is too heavy, and secondly, 
because above 2xl0 19 it has a flat trend, con- 
trary to the measured <ln(A)>, which increases 
with energy. 

Notice that the Auger Collaboration disposes 
today (2009) of a powerful and redundant appa- 
ratus and a number of recorded atmospheric cas- 
cades superior to that of any other experiment. 

The purpose of this work is to calculate the 
chemical composition of the cosmic radiation at 
Earth, not only assuming cosmic-ray sources in 
the disc, but also postulating a component that 
penetrates the disc volume from the exterior 
reaching the Earth. This cosmic-ray component 
is called here extradisc being a particular compo- 
nent of the more vast class of the extragalactic 
cosmic rays. 

Both the extradisc cosmic rays and their flux 
will be denoted by the same symbol I ec i (ed for 
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Figure 2. Energy spectra of individual ions and 
group of ions of Eas-top (red squares) [5] and Kas- 
cade (blue small circles) [10] experiments com- 
pared with the theory of constant spectral indices 
(black curve). Fluxes are normalized at 10 12 eV. 
Ion grouping and flux of the Eas-top experiment 
shown in this figure are examined elsewhere [llj . 



extradisc). Let us anticipate here that extradisc 
component reconciles the results of the theory 
with the Auger data on <ln(A)> shown in fig- 
ure 3, still preserving a quantitative accord with 
the knees, the ankles and the second knee. 

The structure of this paper preassumes the 
existence of the Theory of Constant Spectral 
Indices [1-4] and the companion paper [B] where 
the chemical composition of the cosmic radiation 
is evaluated under the restricted assumption that 
all cosmic-ray sources feeding the particle flux at 
Earth are placed in the disc volume. Section 2 
presents a brief survey of the bases of the calcu- 
lation and a flash of the cultural background of 
the theory. Section 3 is devoted to a survey on 
cosmic-ion energy spectra useful for the normal- 
ization of the theory at 10 12 eV. Section 4 summa- 
rizes the calculation of the chemical composition 
obtained with the galactic sources only. Section 
5 reports a comparison between data and theory 
under the assumption that all cosmic ray sources 
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Figure 3. Comparison of the <ln(A)> derived 
from the Theory of Constant Indices with the 
Tunka [13], Kascade [H] and Auger [T5] data. A 
significant discrepancy between this calculation 
and the Auger data appears above 4 x 10 17 eV. 
The theoretical evaluation of <ln(A)> (blue line) 
assumes that all cosmic-ray sources are placed 
in the disc volume, with no extragalactic com- 
ponent. 



are placed in the disc and not in its exterior. Sec- 
tion 6 reports many details of the calculation of 
the ion abundances of the extradisc component 
and the attenuation of its intensity while pene- 
trating the disc volume from its periphery up to 
the solar system. 

The results of the calculation reported in Sec- 
tion 6 are utilized in Section 7, where the 
<ln(A)> for the disc and extradisc components 
are given. The comparison of the predicted 
<ln(A)> with the data for the two components 
I e d and Id (disc component) are given in Sections 
8 and 9. Section 10 marks the disagreement on 
the chemical composition between the HiRes ex- 
periment and others. Conclusions are in the last 



Table 1 

Parameters of the theory denoted Low Energy 
(LE) and High Energy (HE) ion blends. Fluxes 
are multiplied by E 2 5 and expressed in units of 
(m^sr^s^eV 1 - 5 ). 



Blend. 


LE 




HE 






10 i2 eV 




10 14 eV 






% 


7 


% 


7 


H 


42.4 


2.77 


33.5 


2.67 


He 


26.5 


2.64 


27.3 


2.64 


CNO 


11.9 


2.68 


13.6 


2.65 


Ne-S 


9.2 


2.67 


10.9 


2.63 


Ca(17-20) 


1.2 


2.67 


2.8 


2.63 


Fe(21-28) 


8.7 


2.59 


12.0 


2.62 




Flux 


7 


Flux 


7 


H 


1.15 10 1Y 


2.77 


3.93 10 ib 


2.67 


He 


7.19 10 16 


2.64 


3.20 10 16 


2.64 


CNO 


3.24 10 16 


2.68 


1.60 10 16 


2.65 


Ne-S 


2.50 10 16 


2.67 


1.28 10 16 


2.63 


Ca(17-20) 


3.14 10 15 


2.67 


3.27 10 15 


2.63 


Fe(21-28) 


2.36 10 16 


2.59 


1.41 10 16 


2.62 


Total 


2.71 10 17 


2.70 


1.18 10 17 


2.64 



Section 11. 

2. Bases and foundation of the calculation 

A prerequisite inherent this calculation states 
that cosmic rays generated in the disc volume of 
the Milky Way might migrate into the halo, and 
eventually, into the intergalactic space [6]. The 
disc of Milky Way is regarded as standard, typ- 
ical disc volume. All these cosmic rays escaped 
from the disc volume populating the extradisc 
space, constitute the extradisc component. This 
calculation admits that a fraction of the extradisc 
component existing at the outskirts of the disc 
volume, by displacement predominantly shaped 
by the the magnetic field and gas density, can 
reach the solar cavity (the Earth): its intensity 
will denoted by I e d- 

In order to perform practical calculation it is 
necessary to specify some properties and features 
of the extradisc cosmic rays present outside the 
disc: 
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Figure 4. Energy spectra of 11 individual ions 
measured by balloon and satellite detectors at en- 
ergies below 10 15 eV suggesting and founding the 
postulate of the constant spectral indices ( see ref. 
[16] and references therein). 



Figure 5. Spectral indices of cosmic ions from Hy- 
drogen to Nickel measured by balloon-borne and 
satellite detectors at 10 12 eV (filled red dots) and 
10 14 eV (open blue dots). All indices extracted 
from the spectra around 10 14 eV fall between 2.6 
and 2.7. 



(A) Cosmic rays emitted by the sources placed 
in the disc and escaping from the disc volume un- 
dergo alterations of their abundances as released 
at the sources, due to nuclear interactions with 
interstellar gas. It will be shown that this alter- 
ation depend both on the energy and the type of 
the ion. 

(B) The spectral index of any ion released at 
the sources in the disc volume is constant in the 
energy interval 10 11 -5 x 10 19 eV. The plausibility 
of this fundamental hypothesis is rooted on the 
measurements of the spectral indices collected in 
more than half a century by balloon and satellite 
detectors (see fig. 4, 5 and 6). This hypothesis 



has been converted in 2006 into a postulate called 
Postulate of Constant Spectral Indices pQ. 

(C) Ion abundances as a function of energy of 
the extradisc particles traversing the disc and ar- 
riving at Earth are evaluated assuming the same 
parameters adopted for the disc component Id- 

Though the present calculation ignores the 
mechanism that accelerate cosmic rays in the 
Milky Way, two constraints have to be obeyed by 
the accelerator, whatever it might be, in the range 
10 11 -5 x 10 19 eV : (1) spectral indices generated 
at the sources placed in the disc volume are inde- 
pendent of energy; (2) the cosmic-ray sources are 
uniformly distributed in the disc volume at any 
energy [T2"] . 
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Figure 6. Measured fluxes of individual ions of 
the Wiebel-Sooth Compilation {LE blend) (red 
dots) and HE A Compilation (open blue dots). 
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Figure 7. Theoretical spectra of 6 ions for the LE 
blend and the resulting complete spectrum (thick 
green upper curve). 



It has been demonstrated by full simulation of 
cosmic-ray motion that the displacement of the 
cosmic rays from the sources to any point in the 
Galaxy has a modest effect on the spectral in- 
dices released at the sources. For example, the 
propagation of helium in the disc entails an index 
modification of 0.04 in the energy band 10 12 -10 14 
eV being at the sources 2.65. While for most cal- 
culations this index modification is a negligible 
effect, for the evaluation of the chemical compo- 
sition around the ankle is of quite notable impor- 
tance. 

The interval of application of the Theory of 
Constant Spectral Indices is 10 11 -5x 10 19 eV. The 
upper energy extreme is determined by the real- 
ization that the characteristics of the knee and the 
ankle of the complete spectrum (all-particle spec- 
trum) are quantitatively explained by the theory. 
The theory predicts the exact position where the 
ankle is observed, at (4 — 5) x 10 18 eV and that of 
the knee as well, and many other features. The 
quantitative characteristics of the ankle of the 
complete spectrum derive from those of the an- 
kles of the individual ions. The proton ankle is 
at the lowest energy while the iron ankle is at the 
maximum energy, which turns out to be 5 x 10 19 



eV. 

The theory distinguishes an ankle in the gram- 
mage versus energy and an ankle in the cosmic- 
ray intensity versus energy. The knees are inti- 
mately related by the theory to the ankles of any 
ions, and since there exist an excellent quantita- 
tive agreement with numerous experimental data, 
the maximum energy involved in the knee-ankle 
relationship is the Fe ankle energy. The Fe ankle 
energy in the grammage is at 5 x 10 19 eV, and as a 
consequence, it delimits the range of application 
of this theory. Above the energy of 5 x 10 19 eV the 
results reported in this paper and the companion 
one [5] have to be regarded as numerical exercises. 
In fact, physical phenomena having empirical ev- 
idence operate above 5 x 10 19 eV but they are not 
incorporated in the Theory of Constant Spectral 
Indices as it is presently formulated |l|2|3|4j . 

The results of the calculations depend on the 
geographic site of the sources generating the ex- 
tradisc component. For instance, if the sources 
of the extradisc component would be the ensem- 
ble of the spiral galaxies, due to the similarity of 
the Milky Way Galaxy, only minor modifications 
of the calculations would be necessary to deter- 
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Figure 8. Theoretical spectra of 6 ions for 
the HEA ion blend and the resulting complete 
spectrum (thick upper green curve) normalized at 
10 12 eV with a flux of 3.89 x 10 17 particles/m 2 sr 
s eV 15 . 

mine <ln(A)>. If the sources of the extradisc 
cosmic rays would be placed in the local Super- 
claster of galaxies, the present calculations are 
not sufficiently precise for a number of reasons. If 
the extradisc component are just reentrant parti- 
cles escaped from the disc of the Milky Way, as 
already suggested [2], the calculations are exact. 

Total kinetic energy of particles is used every- 
where in this paper. 

3. Measurements of the spectral indices 
and ion abundances 

Let us preliminarily notice that the interpola- 
tion of the energy spectra of individual ions with 
a constant spectral index in a restricted energy 
band it is an established tradition at energies be- 
low 10 15 eV. An index slightly depending on en- 
ergy would probably conform better to some data, 
in some energy intervals. This refinement, how- 
ever, it is still uncommon, and unnecessary in the 
present work, and consequently omitted here. 

The fundamental hypothesis of the present cal- 
culation incorporated in the Theory of Constant 
Spectral Indices is that all ions at the sources have 



Table 2 

Parameters of the theory at 10 12 eV denoted 
universal and HE A ion blends. Fluxes are 
multiplied by E 2 ^ and expressed in units of 
(m^sr^s^eV 1 - 5 ). 





HE4 




Universal 






% 


7 


% 


7 


H 


42.4 


2.74 


37.2 


2.645 


He 


26.5 


2.72 


26.4 


2.640 


CNO 


11.9 


2.69 


13.9 


2.650 


Ne-S 


9.2 


2.68 


10.1 


2.650 


Ca(17-20) 


1.2 


2.67 


01.5 


2.650 


Fe(21-28) 


8.7 


2.66 


10.8 


2.650 




Flux 




Flux 




H 


1.61E+17 




8.60E+16 




He 


1.01E+16 




6.10E+16 




CNO 


4.54E+16 




3.22E+16 




Ne-S 


3.50E+16 




2.34E+16 




Ca(17-20) 


0.44E+16 




0.34E+16 




Fe(21-28) 


3.31E+16 




2.50E+16 




Total 


3.89E+17 




2.31E+17 





constant spectral indices. 

This simple hypothesis predominantly concurs 
in the prediction of the complete spectrum lead- 
ing to the excellent agreement with experimental 
data as shown, for example, in figure 1. 

Figure 4 reports the energy spectra of 1 1 differ- 
ent ions in the interval 10 10 -10 15 eV. The spec- 
tra form a grid of parallel straight lines in the 
flux units shown. A posteriori the regular forms 
of these spectra with no dips and no spikes jus- 
tify the postulate of Constant Spectral Indices. 
Though the evidence for constant indices is solid 
and irrefutable below 5 x 10 16 eV there exist 
in Nature magnificent and universal phenomena 
that modify indices released at the sources. The 
ensemble of the processes referred to as solar 
modulation occurring in all stars is an example 
of violation of the Postulate of Constant Spectral 
Indices at low energy. Thus, local phenomena ev- 
idently violates the postulate. In order to ascer- 
tain if additional local phenomena in the Galaxy, 
besides solar modulation and the knees, in the in- 
terval 10 11 -5 x 10 19 eV alter the indices, a careful 
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Figure 9. Relative amounts of light, intermediate 
and heavy ions versus energy for the LE (thin 
curves) and HE (thick curves) ion blends. Notice 
the minimum of the light component around 3 x 
10 17 eV in both ion blends, followed by a stable 
increase of the light ion fraction up 4 x 10 18 eV. 

examination of the experimental data is manda- 
tory. 

A compilation of indices of a number of nu- 
clides from Hydrogen to Nickel is shown in figure 
5. This compilation based on measurements of 
the indices mainly between 10 10 -10 15 eV is re- 
ferred to as Wiebel-Sooth compilation \T7\ and 
the corresponding parameters incorporated in the 
theory are called LE ion blend (Low Energy). 
Another compilation of spectral indices collect- 
ing measurements at higher energies |16j between 
10 12 -10 15 eV is also shown in the same figure 5 
(blue open dots)and it is referred to as HE ion 
blend (High Energy) . The spectral indices shown 
in figure 5 accumulate around a common value 
of 2.65 between 10 n -10 15 eV. It is also evident 
from the experimental data on the indices shown 
in figure 5 (red dots) that higher the atomic mass 
harder the index. This correlation between in- 
dices and atomic mass is more pronounced for 
LE than HEA ion blend. 

Notice that the parameters of the Wiebel-Sooth 
Compilation incorporated in the LE ion blend 
based on the interval \QxZ{GeV) up to 10 12 eV 



are three orders of magnitudes from the knee en- 
ergy region (e.g. 10 15 -10 18 eV) while those of the 
HE ion blend are closer. If the postulate of the 
constant spectral indices is exactly described by a 
straight line (in logarithmic scales of intensity and 
energy) either the LE or HE blend would suffice 
to correctly calculate the chemical composition. 
In this circumstance the differences in the chem- 
ical composition at any energy are only due to 
the ion abundances at the normalization energy 
of the theory. If quite small violations of the pos- 
tulate take place, then, the differences originating 
from the two ion blends become significant. 

Recent experimental data between 40 and 400 
GeV indicates that the proton index changes from 
2.77 ± 0.02 to 2.67 ± 0.03 in the interval 4 x 
10 11 4 x 10 14 eV taking into account the data of 
Atic2 [18] and Cream experiments [19]. 

The hardening of the index seems to affect the 
helium spectrum as well. In fact, the index of 
2.70 ± 0.02 measured at low energy below 400 
GeV decreases to the value of 2.62 ± 0.05 above 
400 GeV up to the maximum energy explored to 
day of 4 x 10 14 eV [20]. The alteration of the H 
and He indices in the preknee energy region would 
entail a change in the H/He flux ratio which has 
been proved to be constant at low energy, 10 13 eV 
[2"T] . The H/Hc flux ratio versus energy has been 
examined in another paper [16] where it has been 
shown that the problem of the H/He flux ratio 
above 10 13 eV, is, at the present times, inextri- 
cable. The recent data of the Cream experiment 
[22] do not resolve the problem of the H/Hc flux 
ratio below 10 15 eV due to the large error bars 
of the experiment at high energy, in the decade 
10 14 -10 15 eV. 

Figure 6 reports the ion fluxes for the 2 com- 
pilations defining the parameters of the LE and 
HE blends. The fluxes are almost equal for 
the two blends as shown in fig. 6. The flux 
of the complete spectrum based on the Wiebel- 
Sooth compilation (LE blend) is 2.71 x 10 17 
particles/m 2 s sr eV 15 at 10 12 eV. The HE blend 
at 10 12 is normalized at the same flux which at 
10 14 eV becomes 1.07 x 10 17 particles/m 2 s sr 
eV 1 ' 5 . The complete spectrum has been mea- 
sured in the region 10 12 - 10 15 giving a spectral in- 
dex of 2.74 and an extrapolated intensity at 10 12 



around 2 x 10 17 particles/m 2 s sr eV 1 - 5 [23] . 

In the interval 10 15 -10 17 eV recent flux mea- 
surements of the complete spectrum by Kascade 
give (Sibyll) 7.68 x 10 16 and (QGSjet) 7.48 x 10 16 
particles/m 2 s sr eV 15 at the reference (and ar- 
bitrary) energy of 2 x 10 15 eV. Statistical errors of 
these measurements are about 7 per cent. Note 
that the difference in the complete spectra result- 
ing from QGSjet and Sibyll algorithms is negligi- 
ble (see fig. 1). Differences in flux measurements 
resulting from QGSjet and Sybill algorithms be- 
come quite significant for individual ions amount- 
ing to factors 2-10 for some ions in some energy 
intervals as pointed out elsewhere (see Section 9 
ofref. PI])- 

The Kascade data suggest that the spectral in- 
dices of heavy ions are constant in any preknee 
energy region. Empirically, they are expected to 
be constant on the basis of balloon and satellite 
data on proton and helium spectra collected be- 
low 10 15 eV, provided that the H and He knees 
are regarded as experimentally observed, around 
2 x 10 15 eV and 6.7 x 10 15 eV, respectively. If 
so, heavy ion spectra are expected to behave sim- 
ilarly. 

The heavy ion spectra measured by the Kas- 
cade experiment do have constant spectral indices 
close to 2.65 |24|25j . Any methods of analysis 
(QGSjet, Sibyll or deconvolution method) yield 
indices of about 2.65 or harder but not softer. 
When harder indices are observed a contamina- 
tion between nearby ions is suspected or demon- 
strated (see Section 7 of ref. [H] devoted to this 
argument). With the appropriate ion groupings 
of contaminated adjacent ion samples and the 
correct guide from the theory [1-4], a constant 
common index for heavy ion spectra in the fun- 
damental energy region 10 15 -10 17 emerges [11] , 

Though the present calculation does not de- 
pend on any acceleration mechanism, a theory 

2 In the paper : The Transition from Tortuous to 
Rectilinear Cosmic Ray Trajectories is at the Origin of 
the Knee by A. Codino and F. Plouin [4], the theoret- 
ical Fc spectrum is compared with that measured by the 
Kascade experiment, which has an index of 2.64 ± 0.06. 
In this comparison heavy ions include nuclei from Silicon 
to Nickel. Different groupings result in similar indices but 
always in overabundant fluxes compared to the theory (see 
figs. 16 and 17 of ref. [5] ). 



of cosmic-ion acceleration in pulsar atmospheres 
having constant spectral indices up to 5 x 10 19 eV 
has been conceived in 1969 [26 . 

The findings on the indices of heavy ion spec- 
tra of the Kascade experiment consolidate the 
Postulate of Constant Spectral Indices from 10 15 
eV up to about 10 17 eV (unobserved Fe knee). 
As far as the Postulate of Constant Indices is at 
stake, the outcomes of the Kascade experiment 
on heavy ion spectra in the interval 10 15 -10 17 eV 
have anticipated, probably by more than half a 
century, those plausibly expected by future bal- 
loon and satellite experiments, which identify in- 
dividual ions with superior discrimination power 
but with a very few nuclei. 

The data shown in figure 5 are included in the 
theory with 6 indices and 6 fluxes at the nor- 
malization energy of 10 12 eV. The values of these 
parameters are given in Table 1. 

In the same Table 1 are given the parameters 
of the HE blend tuned to the data at higher 
energy. In the following, in order to appreciate 
differences and similarities with previous calcula- 
tions, the parameters of 2 additional ion blends 
referred to as universal and HEA blends are given 
in Table 2. In the universal ion blend all in- 
dices are equal with a value of 2.65 except for 
helium set at 2.64 which takes into account the 
recent trend of the measurements suggesting an 
index harder than that of other ions. The HEA 
blend is a variant of the HE blend where the 
dominant nuclei, e.g. protons and helium, H and 
He have the classical value of 2.74 [53] and 2.72 
[2"7] believed for many years to be standard, reli- 
able measurements. The parameters of the HEA 
blend are used to determine the chemical compo- 
sition around the ankle (see fig. 6). The HEA 
blend is normalized at 10 12 eV with a flux of 
3.79 x 10 17 particles/m 2 s sr eV 1 - 5 . This value 
is inspired from the Kascade flux measurement 
which is 7.58 x 10 17 particles/m 2 s sr eV 1 ' 5 at 
2 x 10 15 eV (mean value of QGsjet and Sibyll algo- 
rithms) [7]. The HEA flux normalization at 10 12 
eV is larger than that of the LE blend (Wiebel- 
Sooth flux compilation) by a factor 1.4. 

Because of the ion propagation in the disc vol- 
ume the indices of the original spectra released 
by the sources modify (see fig. 6 of ref. 3). 



It may be useful to subdivide the index of any 
ion spectrum measured at the solar cavity in two 
parts: the propagation index and the accelera- 
tor index at the sources. The propagation indices 
due to ion displacement in the disc volume for 
the 6 ions H, He, N, Si, Ca and Fe are, respec- 
tively, 0.060, 0.055, 0.045, 0.045, 0.040 and 0.040. 
Taking into account the accelerator indices of the 
universal blend (see Table 2) the global index be- 
come: 2.585, 2.585, 2.605, 2.605, 2.610 and 2.610 
for the same ion sequence just mentioned. 

4. Chemical composition of the cosmic ra- 
diation generated by cosmic-ray sources 
placed in the disc 

The energy spectra of 6 ions derived from the 
Theory of Constant Spectral Indices with the LE 
and HEA blends are, respectively, in figure 7 and 
8. The complete spectrum which is the sum of the 
partial spectra of individual ions is also shown in 
the same figures 7 and 8 for both blends result- 
ing in a flux of 8.2 x 10 16 particles/m 2 s sr eV 1 ' 5 
for the HE4 ion blend at 10 15 eV. The ion abun- 
dances derived from the energy spectra reported 
in figures 7 and 8 are shown in figure 9 for the 
interval 10 12 -5 x 10 19 eV. The light ion fraction 
(H+He) decreases with a gentle slope going from 
68 per cent at 10 12 eV to 55 per cent at 1.5 x 10 15 
eV (HE4 ion blend) . Above this energy the light 
ion fraction descends abruptly reaching a mini- 
mum of 11 per cent at the energy of 4 x 10 17 eV, 
quite evident in figure 9. A similar behaviour is 
exhibited by the LE blend but with higher light 
ion fractions. 

The slight increase in the (Ca+Fe) fraction in 
the interval 10 12 -4 x 10 15 eV for the HEA blend 
is caused by the tiny difference in the Ca and 
Fe slopes (2.66 for Fe and 2.650 for Ca, see Ta- 
ble 2) with respect to the light ion slopes (2.74 
for H and 2.72 for He). The strong increase of 
the (Ca+Fe) fraction between 4 x 10 15 eV and 
7 x 10 17 eV is due not only to the difference in 
the indices(2.70 light ions against 2.60 heavy ions) 
but predominantly to the high fall of the light ion 
flux (light ion knees) above the nominal knee of 
the complete spectrum at 3 x 10 15 eV. The fall is 
quite notable for protons and helium in the inter- 



val 2 x 10 15 -2 x 10 17 eV but insignificant for heavy 
ions (Ca+Fe) at least up to 4 x 10 17 eV. The pro- 
cesses causing this abrupt fall of intensity (knees) 
have been analyzed and described elsewhere [3]. 

5. The <ln(A)> of cosmic-ray disc compo- 
nent compared with the data 

It is instructive to compare the theory with 
the experimental data on <ln{A)> under the re- 
stricted assumption that all cosmic-ray sources 
are in the disc volume with no extragalactic com- 
ponent. 
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Figure 10. Mean logarithmic mass of the cosmic 
radiation resulting from the Theory of Constant 
Indices in the energy interval 10 n -5 x 10 19 eV 
for three ion blends. In these 3 evaluations all 
cosmic-ray sources are placed in the disc volume. 

The ion spectra shown in figure 7 or 8 directly 
determine the chemical composition of the cosmic 
radiation in terms of </n(j4)>. The <ln(A)> at 
energy E is the weighted average of the logarith- 
mic atomic masses of the cosmic rays: 

<ln(A)>= £ f A ln(A) 

where A is the atomic mass of the nucleus and 
fA its fraction in the cosmic-ray flux. 

In figure 10 are shown 3 profiles of <ln{A)> re- 




Figure 11. Comparison of the theoretical 
<^n(yl)> for the LE and HE4 blends with that 
measured by Kascade (black triangles) [28] using 
the deconvolution method and that measured by 
Eas-top (red dots) [29] . 

suiting from 3 different ion blends: LE, HEA and 
the universal blend. The differences in the LE 
and HEA blends are quite negligible in the entire 
interval 10 12 -5 x 10 19 . There are significant sim- 
ilarities but also some differences. The 3 profiles 
are similar in the band 2 x 10 15 -2 x 10 17 eV, they 
all increase by about 2 units of <ln{A)> and, in 
the range 5 x 10 17 -5 x 10 18 , they all descend by 
about one unit. Differently, the <ln(A)> profile 
of the universal blend is flat from the initial en- 
ergy at 10 12 up to 2 x 10 15 eV, then it increases 
reaching a maximum of 3.05 at 2.97 x 10 17 eV, 
and a more pronounced descent at high energy, 
and then, above 5 x 10 18 eV, a modest decreasing 
trend dominates. The form of the <ln{A)> pro- 
file of the HEA blend increases up to a maximum 
value of 3.33 at 3.94 x 10 17 eV, then it decreases 
rapidly up to (4—5) x 10 18 eV leveling off at higher 
energies around 2.8. As explained elsewhere (see 
fig. 4 of ref. pQ) the spectra of heavy ions (Ca 
and Fe) attain the low asymptotic plateaux (see 
fig. 7 and 8) at higher energies than the light 
ions. Low asymptotic plateaux reflect also the 
quasi rectilinear propagation in the Milky Way. 

Figures 11 and 3 show data on <ln(A)> of the 
Kascade experiment which separates all nuclei in 



Figure 12. Atmospheric depths of shower maxima 
of protons and Fe nuclei versus energy computed 
by QGSjet and Sibyll hadronic codes as evaluated 
by Heck [30]. 

5 groups along with the corresponding theoret- 
ical curves (LE and HEA). The Eas-top data 
on <ln(A)> are also shown in figure 11. Though 
there are discrepancies in the ion spectra in figure 
3 obtained by QGSjet and Sibyll algorithms of the 
Kascade experiment [TT], the resulting <ln(A)> 
profiles show the correct general trends also ob- 
served by other experiments. By correct trend is 
meant an average increase of 1.5 units of <ln(A)> 
in the interval 10 15 -10 17 eV from an initial value 
of 1.8 at 10 15 eV dictated by the extrapolation 
of balloon and satellite data at lower energies. 
The forms of the <ln(A)> profiles in figure 11 
obtained by the method of deconvolution [?] do 
not disagree with those of QGSjet and Sibyll al- 
gorithms shown in figure 3, though the <ln(A)> 
in some energy bands have notable differences for 
the 3 methods. The agreement of the theoreti- 
cal <ln(A)> profiles for the LE and HEA blends 
with the Kascade data has been described in de- 
tail elsewhere [BJ. 

At energies larger than 10 17 eV the determi- 
nation of chemical composition of the cosmic ra- 
diation becomes more involved because ions are 
not resolved individually nor in restricted groups. 
In order to determine the <ln(A)> from the ob- 
servables of the atmospheric showers recorded by 
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Figure 13. Atmospheric depths of shower maxima 
X m ax versus energy measured by Fly's Eye (blue 
triangles) [31] and Yakutsk (red stars) [32] exper- 
iments compared with the corresponding theoret- 
ical profile (red line). 



detectors, accurate and detailed models of the 
nucleus-air interactions are required. Tradition- 
ally, the depth in g/cm 2 of the maximum of the 
atmospheric showers (i.e. number of particles ver- 
sus cascade axis) is denoted X max or equivalently 
elongation. 

A number of methods have been devised to 
measure X max : (1) muon number and muon den- 
sity detected at ground; (2) the mean width of the 
X m ax distribution denoted in short a(X max ); (3) 
Cherenkov light generated in the cascades; (4) flu- 
orescence light generated in the cascades; (5) cur- 
vature of the cascade front measured by time-of- 
flight techniques. Apparata exploiting Cherenkov 
and fluorescence light can be equipped with ad- 
ditional detectors. 

Figure 12 gives the X max versus energy evalu- 
ated by Heck [30] for two models of nuclear in- 
teractions denoted in short QGSjet and Sibyll. 
It results that the depth of the maximum of the 
atmospheric showers induced by primary cosmic 
proton goes from 500 g/cm 2 at 10 14 eV to 850 
g/cm 2 at 10 18 eV. In the same energy interval 
the Fe elongation goes from 350 g/cm 2 at 10 14 
eV to 750 g/cm 2 being the total thickness of the 
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Figure 14. Atmospheric depths of shower maxima 
of protons and Fe nuclei versus energy measured 
by Auger, Space2/Vulcan [53] and Casa-Blanca 
[34] experiments compared with the <ln(A)> 
profile derived from the theory. 

air about 1000 g/cm 2 . 

Above 10 14 eV the gap between the proton 
X max profiles of Sibyll and QGSjet regularly en- 
larges reaching 40 g/cm 2 at 10 20 eV (see fig. 12). 

According to some simplifications the elonga- 
tion is given by: 



X,, 



D((lnE/E ) -ln(A))+C 



where X rnax is the atmospheric depth in g/cm 2 , 
for primary particles of energy E, Eq is a refer- 
ence energy, <ln(A)> the chemical composition 
and C and D two appropriate functions. Let us 
note that in the present calculation, unlike others, 
C and D vary with energy. In the following, in 
the comparison with the experimental data, both 
X max and <ln(A)> are maintained since the two 
variables are affected by small differences. The 
function X^ ax (E,Eo) is the elongation for the 
nucleus A with the energy E normalized at the 
energy E . 

The function X max (E,E ) obtained by cascade 
simulations depends on the hadronic models used 
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Figure 15. Chemical composition in terms of 
<ln(A)> extracted by the X max measured by 
Casa-Blanca |34 and Auger experiments. The 
<ln(A)> profiles of each experiment curiously ex- 
hibit prominent depressions at 4 x 10 15 eV and 
(2 — 3) x 10 18 eV. The Casa-Blanca depression is 
probably due to an instrumental effect as argued 
in the text. 

to describe the interactions of cosmic nuclei with 
the air. Once the experimental value of the 
elongation X^ x has been determined, the 
value of <ln(A)> is obtained in the superposi- 
tion model by the equation: 



<ln(A)>- 

max) 



Zn(56)x (X% ax -X%P x )/(X% a 



being X^ ax and X^% x the elongations sim- 
ulated by Monte Carlo for protons and Fe nu- 
clei, respectively. Using the above equation, the 
<ln(A)> derived from the theory can be con- 
verted into X max . The result is shown in figure 
12 with the blue and red curves which are the av- 
erage value of X^ax an< ^ ^mlx profiles obtained 
by QGSjet and Sibyll codes. 

The Yakutsk [32] and Fly's Eye [3T] data on 
<ln(A)> are shown in figure 13 and they rep- 



resent an example of accord between theory and 
data. Note that at energies above 10 18 eV the the- 
oretical X max profile has a substantial gap with 
the data. The accord below 10 18 eV and the dis- 
crepancy at high energy have been discussed in 
the companion paper [6]. 

Figure 14 reports examples of disagreement. 
The X max measured by Auger [15] above 10 17 
eV exhibits a large discrepancy between data and 
theory in the huge interval 4 x 10 17 - 10 19 eV. Sub- 
sequent elaborations of the Auger data samples 
[?] mark the stability of the X max except for the 
last data point which decreases by about 5 g/cm 2 . 
The Casa-blanca and Space2 /Vulcan data points 
in figure 14 have opposite deviations with respect 
to the theory, in different energy intervals. 

The measurements of the atmospheric depths 
of the cosmic radiation made by Auger has been 
converted into <ln{A)> and shown in figure 15. 
According to these measurements cosmic rays 
around 4 x 10 18 eV mainly consist of 37 per cent 
of heavy ions. This figure comes from the parti- 
tion of the cosmic nuclei in two groups, light and 
heavy (see subsequent Section 8). 

Besides the Auger experiment, a notable struc- 
ture in the <ln(A)> profile, shown in figure 15, 
has been also observed by the Casa-Blanca ex- 
periment in the interval 10 14 - 4 x 10 16 eV. The 
depression of <ln(A)> measured by Casa-Blanca 
in the energy band (3 — 9) x 10 15 eV would im- 
ply the disappearance of heavy elements. This 
depression corresponds to the decreasing differ- 
ence between the measured X max profile and the 
theoretical X" ax in the interval 2 x 10 14 -3 x 10 15 
eV as shown in figure 14. As a consequence, the 
depression of the Casa-Blanca experiment shown 
in figure 15 is not a spurious effect in the conver- 
sion of X max into <ln(A)>. This circumstance 
favours an instrumental effect and not a physical 
cause for the origin of the depression. The con- 
clusion is corroborated by the measured by the 
Kascade Collaboration with 3 methods (fig. 3 
and 11) and others experiments which, contrary 



to Casa-Blanca data, register the disappearance 
of light elements. 







6. Ion abundances outside the disc volume 



In order to determine the ion abundances in 
the space outside the disc volume it is necessary 
to calculate the probability of escaping from the 
disc, Pe, for any individual ions at all energies. 
At very high energy, as ions propagate almost rec- 
tilinearly, this calculation is extremely simple. 
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Figure 16. Probability of escaping from the disc 
Pe as a function of the energy for particles hav- 
ing cosmic-ray sources distributed uniformly in 
the disc volume [12] ■ The probability Pe is nor- 
malized to 2 x 10 5 particles. 

Let E d A be the energy above which ions of mass 
A propagate almost rectilinearly in the disc. In 
the rectilinear propagation regime ions encounter 
a minimum amount of matter as evaluated else- 
where [3] . 

The amount of matter experienced by ions 
while propagating in the disc is referred to as 
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Figure 17. Number of nuclear collisions versus 
energy occurring in the disc volume for nuclei (H, 
He and Fe) having all the initial points of the 
trajectories placed in the solar cavity. 

grammage , g. It is given by: g=mnL, where m is 
the mean atomic mass in the disc, n is the aver- 
age number of atoms per cm 3 in the interstellar 
space and L the trajectory length. The grammage 
depends on nuclear cross sections via trajectory 
length, L. Since for E > E\ the grammage is 
only 0.006 g/cm 2 , the number of nuclear collisions 
taking place in the disc is negligible, and conse- 
quently, ion abundances released by the cosmic- 
ray sources are almost unmodified when they are 
observed in the solar cavity. 

Unlike the rectilinear propagation regime in the 
energy interval where E < E A , between the knee 
and the ankle energy region, e. g. 3xl0 15 -5xl0 18 
eV, the galactic magnetic field affects ion mo- 
tion, bending and inverting efficiently ion direc- 
tions. This effect of the magnetic field for each ion 
makes the average grammage encountered by cos- 
mic rays dependent on particle momentum. For 
E < E A , the calculation of Pe and the proba- 
bility for cosmic ions of entering (or re-entering) 
from the disc boundary to the solar cavity, Pr, is 
more involved and difficult. In the following the 



functions Pr are determined. 



The number of particles escaping from the disc 
Np versus energy for 6 ions or group of ions is 
given in figure 16. The quantity Np is normalized 
at the same number (2 x 10 5 ) of emitted particles 
(becoming the probability Pe). All curves in fig- 
ure 16 have a minimum at a particular energy 
E" lm (e is for escape). The characteristic ener- 
gies where the minima of the 6 ions H, He, N, 
Si, Ca and Fe occur are, respectively : 1.1 x 10 15 
eV, 2.6 x 10 15 eV, 1.2 x 10 16 eV, 2.4 x 10 16 eV, 
3.9 x 10 16 eV, and 4.5 x 10 16 eV. The decreas- 
ing segment of Np from 10 12 eV up to E™ m is 
due to the increasing nuclear cross sections with 
energy while the rising segment above E™ %n is 
predominantly caused by the vanishing efficiency 
of the galactic magnetic field with rising energy 
to retain particles for long times. 

Escaped particles populating the extradisc 
space may re-enter the disc in the scheme delin- 
eated in Section 2, in a variety of circumstances. 
Let be Pr the mean probability for ions of energy 
E of reaching the Earth from the disc boundary. 
A cosmic-ray trajectory consists of an initial point 
(the source) and a termination point, which is the 
location where the nuclear collision occurs within 
the disc volume. The calculation of Pr exploits 
the approximate symmetry of the cosmic-ray tra- 
jectories contained in the disc volume between the 
initial point Xj and the final point Xp of the tra- 
jectory. Samples of 2 x 10 5 particle are injected 
from the solar system position and the trajecto- 
ries reconstructed in the disc volume. Let Ni be 
the number of ions stopping in the disc volume 
by nuclear collisions and Nt the total number of 
trajectories (or injected particles). The function 
Nj is of fundamental importance in the present 
calculation since it incorporates the effect of the 
magnetic field, the dimension of the disc, the vari- 
ation of the nuclear cross sections with the energy 
and the position of the solar system in the Galaxy. 
The calculation is finally accomplished by taking 
advantage of the inversion of the particle trajecto- 
ries, the interchange of Xp with Xj 7 as explained 
in detail elsewhere [35] ) 
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Figure 18. Probability of reaching the solar sys- 
tem for cosmic nuclei having all the initial points 
of the trajectories placed at the disc boundary. 

Figure 17 reports the functions Ni for 3 rep- 
resentative ions H, He and Fe taken as reference 
examples. The Ni profile of any ion has a high 
plateau, a quite small rise controlled by nuclear 
cross sections, a rapid and high descend, a min- 
imum, and finally, a modest increase up to the 
highest energies. For example, the number of He 
interactions in the disc, in the interval 5 x 10 11 - 
10 15 eV, increases by 3.5 per cent. At 10 12 eV 
the number of He particles interacting in the disc 
volume is 173000 out of 200000 injected from the 
Earth (hence, N T =2 x 10 5 and AT Z =1.73 x 10 5 ). 
At 5 x 10 18 eV where the minimum of Ni for He 
nuclei occurs, the number of He interactions re- 
duces to 2000 (hence iV/=2000). Therefore, the 
probability for He particles of reaching the solar 
cavity, having the initial points of the trajecto- 
ries at the disc boundary, is : Pr = 1 - (iVj/iVy). 
With the figures given above, the probability Pr 
is 0.11 at 10 12 eV and 0.99 at 10 18 eV. By this 
procedure the 6 probability functions versus en- 
ergy Pr are obtained and shown in figure 18. The 
number of nuclear collisions taking place in the 
disc volume shown in figure 17 (see figure 6 of 
ref. [T]) is converted in terms of probability Pr 
in figure 18. The normalization of the functions 



in figure 17 is explained in detail on the investi- 
gation on the origin of the knee (Section 4 of ref. 
[3]). Thus, the functions in figure 17 are closer to 
the simulation of the cosmic-ray trajectories and 
they facilitate the comprehension of the processes 
causing the ion knees [3], while those in figure 18 
more specifically relate to the present study. 
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Figure 19. Particle fraction, f, versus energy of 
the extradisc component reaching the solar cav- 
ity from the disc boundaries (Halo of the Milky 
Way). The quantity / is defined as f—N e T where 
N e is the number of particles escaped from the 
disc and T the transmission probability (see fig. 
18) from the disc boundaries to the solar cavity 
at a specified energy. 

Figure 20 is an illustration of how ion abun- 
dances of the extradisc component are affected 
by small variations of the spectral indices at the 
sources for cosmic ions migrating from the disc 
into the extradisc space. The abundances of H, 
He and Fe ions in the extradisc volume at 10 19 
eV are taken with the arbitrary ratios of 1.5, 1 
and 0.33 (see fig. 20). Note that the three black 
curves (H, He and Fe) in fig. 20 are the same 
reported in figure 16 except for a different nor- 
malization. The Helium index jjje is taken as 
a reference value (any value, including 7jj e = 0) 



.8 




Energy (eV) 



Figure 20. Small differences of ±0.1 % in the 
spectral indices at 10 12 eV will reverberate large 
alterations of the ion abundances around the an- 
kle energy region. The grid of red and blue 
curves, specifying H and Fe ion abundances ver- 
sus energy, result from the change, by arbitrary 
steps of 0.02 of the H and Fe indices, with respect 
to the He index taken as reference value (see text 
for details). 

while the Fe and H indices are varied by small 
amounts of 0.02, 0.04, etc. with respect to 7# e 
, in order to map how ion abundances vary with 
energy. The result is displayed in figure 20 by 
the grid of red and blue curves in the interval 
10 12 -5 x 10 19 eV. For example, with a change of 
0.08 in the index, the H/Fe abundance ratio of 
19.9 (arbitrary value) at 10 12 eV becomes 0.23 at 
10 19 eV, a factor 86.5 minor. Such a figure rep- 
resents a large variation of the chemical compo- 
sition. Note that at very high energy, above 10 19 
eV, the abundance ratios outside the disc volume 
are the same existing at the cosmic-ion sources in 
the disc, since particle displacement is unaffected 
by nuclear collisions, due to the small grammage 
traversed. 

Measurements of the spectral indices in the en- 
ergy region 10 11 -10 15 eV acquire a fundamental 




Figure 21. Energy spectra of the extradisc ions 
reaching the solar cavity emitted by sources char- 
acterized by the LE blend. The black star defines 
the flux normalization. 



Figure 22. Energy spectra of the extradisc ions 
reaching the solar cavity emitted by sources char- 
acterized by the HEA ion blend. The black star 
signals the flux normalization. 



importance because their values determine the 
chemical composition at very high energy. The 
postulate of Constant Spectral Indices, incorpo- 
rating the balloon and satellite data below 10 15 
eV, relates the chemical composition at low en- 
ergy to that at very high energy by logical neces- 
sity avoiding any ad hoc mechanisms, at small, 
arbitrary energy bands, to explain the chemical 
composition of the cosmic radiation. 



7. The Chemical composition of the ex- 
tradisc component 

The chemical composition of the extradisc com- 
ponent is radically different from that of the disc 
component. Heavy ions, escaping from the disc 
volume and penetrating through the disc of the 
Milky Way from its exterior, pay twice a trib- 
ute to larger nuclear interaction cross sections 
with respect to light ions. Therefore, the heavy- 
to-light ion flux ratio inevitably augments above 
10 15 eV up to the ankle energy region. 



Figures 21 and 22 report the energy spectra 
of the extradisc ions reaching the solar cavity 
for the two LE and HE4 blends. Firstly, note 
that a simplified situation takes place above the 
maximum of the Fe energy spectrum located at 
2xl0 19 eV. The ion abundances above this en- 
ergy coincide with those at 10 12 eV reported in 
Table 1. Above 2xl0 19 eV the grammage is so 
small that nuclear interactions in the interstellar 
medium do not alter the chemical composition 
of the extradisc component existing at the out- 
skirts of the disc volume. The complete spectrum 
in figure 21 and that in figure 22 are normalized 
at the arbitrary energy of 10 19 eV to the galac- 
tic flux of 5.35xl0 14 (LE) and3.91xl0 14 (HE4) 
particles/m 2 sr s eV 15 . 

Figure 23 shows ion fractions versus energy of 
the extradisc component for 3 groups of ions: 
light, intermediate and heavy being, respectively, 
(H+He), (CNO + Ne-S) and (Ca + Fe). This ion 
partition is adopted in some experiments. Let 
us now analyze the ion fraction profiles starting 
from the extreme high energy where the situation 




Figure 23. Relative amounts of light, interme- 
diate and heavy ions versus energy for the LE 
(thin curves) and HEA (thick curves) ion blends 
generated by the extradisc component. 

is simpler. 

Just above 2 x 10 19 eV ion abundances of the 
extradisc component reaching the solar system 
are equal to those of the cosmic-ray component 
present in the parent galaxy, because particle dis- 
placement through the disc (escaping and enter- 
ing) do not alter ion abundances. Around 2 x 10 19 
eV terminates the increase of 2 orders of magni- 
tude of the (Ca+Fe) ion fraction. The (Ca+Fe) 
fraction above 5 x 10 18 eV does not surpass 17 
per cent, even with an increase of a factor 66 from 
the low level flat fraction, below 5 x 10 16 eV. The 
intermediate ions, i.e. the sum of CNO and (Ne- 
S) group of nuclei, has a similar pattern with a 
step, between the low and high level, of a factor 
25, which constitutes 31 per cent of the extradisc 
component above 5 x 10 18 eV. 

From the profiles of the ion fractions shown in 
figure 23 or from the more detailed ion spectra 
shown in figures 21 and 22, the <ln(A)> of the 
extradisc component is directly calculated; it is 
shown in figure 24 (pink curve, HE A blend). The 



Figure 24. Chemical composition in terms of 
<ln(A)> versus energy of the disc (blue curve) 
and extradisc component (pink curve) for the 
HEA ion blend derived from the Theory of 
Constant Spectral Indices. 

same figure 24 shows, for visual comparison, the 
corresponding <ln{A)> resulting from the disk 
component (blue curve, blend HEA), which has 
a quite different pattern. The <ln(A)> for the 
LE blend of I e d has a quite similar pattern (fig. 
10). The two <ln(A)> profiles of the disc and the 
extradisc components in figure 24 are normalized 
at the same cosmic-ray intensity at Earth. 

The <ln(A)> of the extradisc component ex- 
hibits a constant value below 5 x 10 16 eV, fol- 
lowed by a rising trend in the interval 5 x 10 16 - 
2 x 10 19 eV with an approximate constant slope. 
Above 2 x 10 19 eV the <ln{A)> stabilizes to 
the constant value of 2.1. The physical phenom- 
ena shaping this silhouette have already been dis- 
cussed. A simplified, instant, qualitative compre- 
hension might be depicted in terms of grammage 
encountered by the cosmic rays during the dis- 
placement in the Galaxy. When the grammage is 
large (see figure 16 ref. [3]) heavy nuclei are de- 
stroyed, when the grammage is minimum (0.006 
g/cm 2 ) heavy nuclei travel undisturbed like light 
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Figure 25. Chemical composition in terms of 
<ln(A)> versus energy for various values of the 
extradisc-to-disc flux ratios, r, for the HE4 ion 
blend. 

nuclei and the persists unchanged regardless of 
the energy. In the intermediate region 5 x 10 16 - 
2x 10 19 nuclear cross sections and geometrical fac- 
tors compete and interplay, forging the chemical 
composition with the stable, rising trend shown 
in figure 24. 

When the disc and extradisc components 
are combined together for different values of 
f=I e d/Id, it results the surprising <ln(A)> pro- 
files shown in figure 25. For r=l the flux in the so- 
lar cavity is 1.8 x 10 14 particles/ (to 2 sr s eV 15 ) at 
10 19 eV for both disc and extradisc components. 
The surprising aspect of the <ln(A)> profiles in 
figure 25 resides in the fact that they thoroughly 
resemble to the profile of the experimental data 
on <ln(A)> extracted from the measurements of 
X m ax of the Auger experiment. 

Flux measurements of the cosmic radiation at 
Earth imperatively constraint r, the maximum 
value of the extradisc component, as evident in 
figure 26. The sum of the fluxes of the disc and 
extradisc components is shown in figure 26 along 



Figure 26. Energy spectra of the cos- 

mic radiation with different values of the 
extradisc-to-disc flux ratio, r, derived from the 
Theory of Constant Indices framed between 
the Haverah Park and Auger flux measurements. 



with flux data of two experiments, taken as an 
example. Magnifying the extradisc component to 
match the Auger data on <ln(A)> would entail 
a conflict with the observed fluxes. 

Notice that the theoretical spectra in figure 26 
intrinsically exhibit the adequate change in the 
slope around 4 x 10 18 eV, from 3.2 to 2.7, which 
constitutes the ankle. Similarly, around 10 15 eV, 
the slope of the complete spectrum changes from 
2.7 to about 3.0 which is a characteristic features 
of the knee. 

With the assumptions adopted in this paper 
any types of extradisc component generate a 
rather light chemical composition below 10 17 eV 
due to the dominant role of nuclear cross sections 
suffered by heavy ions while propagating from pe- 
riphery to disc core. 
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Figure 27. Comparison of the theoretical profiles 
of <ln{A)> (red and pink curves) with the em- 
pirical one extracted from the X max measured by 
the Auger experiment (red diamonds and dots). 
The extradisc-to-disc flux ratio at 10 19 is 2 (blue 
line) and 10 (pink line). 

8. Comparison between theoretical and 
measured <ln(A)> with the disc and ex- 
tradisc components 



In the following the experimental data on 
<ln(A)> are compared with the profile derived 
from the theory in the condition r=2. Note that 
the <ln(A)> in the the condition r=l or r=3 
would describe the experimental data as well, 
due to the large systematic uncertainties inherent 
to the hadronic models describing atmospheric 
showers. The value r=2 is regarded here as a 
adequate compromise between cosmic-ray inten- 
sity (see fig. 26) and the depth of the minimum 
of the <ln(A)> profile (see fig. 25). 

In all the subsequent figures 27, 28, 29, 30, 31, 
32, and 33, where high energy data on <ln(A)> 
above 10 17 eV are examined and compared with 
the theory, the <ln(A)> determined by the Kas- 
cade experiment in the range 10 15 -10 17 eV is dis- 



A 
< 

V 




V 2 io 13 io 14 io 15 io 16 io 17 io 18 10 1 



Energy (eV) 



9 io 20 



Figure 28. Comparison of the theoretical profile 
of <^n(j4)> (red curve) with the empirical one 
extracted from the data of Fly's Eye experiment 
(blue triangles). 

played as adequate, reference measurement. 

From the Kascade Collaboration three differ- 
ent procedures to determine <ln(A)> are avail- 
able called here Sybill algorithm, QGSjet algo- 
rithm and deconvolution method. Data of the 
first two methods are in figure 3 and 27 while 
those of the deconvolution method, which covers 
the largest range, i.e. 9 x 10 14 -1.4 x 10 17 eV, are 
in figure 28. The resulting <ln(A)> profiles from 
the three methods of measurement by Kascade 
are simultaneously shown elsewhere (see figure 24 
of ref. [11]) exhibiting a coeherent silouhette in 
the entire interval 10 15 -10 17 eV. 

Recent measurements of X max of the Auger 
Collaboration converted into <ln(A)> by the ref- 
erence profiles for H and Fe shown in figure 14 are 
given in figure 27. Two predicted <in(^4)> pro- 
files according to the Theory of Constant Indices 
with r—2 and r=10 are shown in fig. 27. In 
the condition r=2 the Auger data in the interval 
4 x 10 17 -10 19 eV are still below the theoretical 




Figure 29. Comparison of the theoretical profile 
of <ln(A)> (red curve) with the empirical one 
extracted from the data of HiRes experiment be- 
fore (blue data points)and after (red data points) 
acceptance corrections. 

profile. 

The <ln{A)> extracted from the X max mea- 
sured by the Fly's Eye experiment [31] is shown 
in figure 28. The increasing trend of <ln(A)> in 
the interval 10 17 -6 x 10 17 eV contrasts with the 
opposite trend of HiRes and Agasa-Akeno exper- 
iments (see figs. 29 and 32); it also disagrees with 
a flat silouhette observed by Yakutsk in the same 
energy interval, and with a value of 3.5 measured 
by Kascade. The behavior of </n(j4)> in figure 

28 is reminescent of that observed by the Hav- 
erah Park experiment (see fig. 31) in the same 
energy interval, where a lack of uniformity in the 
detector acceptance is suspected. 

The profile of <ln(A)> extracted from the 
X ma x measured by the HiRes experiment is 
shown in figure 29. The HiRes data in figure 

29 include HiRes Prototype [36] , HiRes Stereo 
[37] and a data revision which takes into account 
acceptance correction and detector performance 



Figure 30. Comparison of the theoretical profile 
of <ln{A)> (red curve) with the empirical one ex- 
tracted from the data of the Yakutsk experiment 
(blue small circles). 

38 . Around 4 x 10 18 eV data reported in fig- 
ure 29 have a mean value of <ln(A)> of about 
1. This tiny value would correspond to the min- 
imum observed by Auger at (2 — 4) x 10 18 eV 
as far as a vague increase of <ln(A)> would de- 
lineate in the band 5 x 10 18 -4 x 10 19 eV. The 
chemical composition extracted from the X max 
measured by HiRes is lighter than that resulting 
from the highest permissible values of the param- 
eter r (dominance of extradisc component) in the 
interval 10 17 -10 19 eV. It is also lighter, by about 
one unit of <ln(A)>, than that of Auger, Fly s' 
Eye, Haverah Park and Agasa experiments (see 
figure 4 of ref. [39].) 

The profile of <ln(A)> extracted from the 
X max measured by the Yakutsk Collaboration 
[32] is shown in figure 30. There is a global, re- 
markable agreement with the theory in the huge 
interval 10 15 -10 19 eV characterized by a bell- 
shaped profile with a maximum of <ln(A)> of 3 
around the energy of 2 x 10 17 eV. Above 5 x 10 18 
eV data would suggest the lightning of the chem- 




Figure 31. Comparison of the theoretical profile 
of <ln(A)> (red curve) with the empirical one 
extracted from the data of the Haverah Park ex- 
periment (blue triangles). 

ical composition. Notice, by contrast, that the 
chemical composition of measured by Auger in 
the interval 5 x 10 18 -4 x 10 19 eV would indicate 
the opposite trend: the <ln(A)> becomes heav- 
ier according to the data shown in figure 27. 



9. Comparison between theoretical and 
measured <ln(A)> in Haverah Park, 
Agasa and Volcano Ranch 

The partition of the cosmic nuclei in two 
groups, light and heavy, with fractions F p and 
Fp e respectively, characterizes the chemical com- 
position extracted from Haverah Park , Agasa 
and Volcano Ranch data samples. This bi- modal 
data analysis, though appropriate for some de- 
tectors, has a clear bias since it gives intrinsically 
higher <ln(A)> than realistic ion groupings as 
demonstrated below. 

The chemical composition measured by the 



Figure 32. Comparison of the theoretical profile 
of <ln(A)> (red curve) with the empirical one ex- 
tracted from the data of Agasa experiment (blue 
and pink data points). 

Haverah Park experiment is reported in figure 
31. Ruling out the data points at the extreme 
energies, below 2 x 10 17 eV and above 10 18 eV, 
where detector acceptance for atmospheric show- 
ers might have been deteriorated (See Section 3 
ref. [10]). the agreement with the theoretical 
<ln(A)> profile in the condition r—2 is excellent. 
For protons and Fe nuclei the differences in the 
theoretical X max profiles shown in figure 12 for 
Sibyll and QGSjet codes arc small in the interval 
10 17 -10 18 eV, and so the agreement between the- 
ory and Haverah Park data is quite stable. It has 
been remarked [41] the critical role of hadronic 
models in the partition of all cosmic-ray nuclei in 
two classes, light (F p ) and heavy (Fp e ). Adopting 
the code QGSjet98 it results a proton fraction F p 
of .34 while the code QGSjetOl gives F Fe = 0.48 

m- 

Data samples collected by the Agasa Collabo- 
ration have been revisited in 2003 [42] in order to 
determine the X max of the cosmic radiation in the 
range 10 17 -10 20 eV. The chemical composition re- 




Figure 33. Comparison of the theoretical profile 
of <ln{A)> (red curve) with the empirical one 
extracted from the data of Volcano Ranch experi- 
ment (blue cross in the rectangle). The data point 
corrected according to the Theory of Constant 
Spectral Indices is represented by a pink cross in 
the rectangle. 

suiting from this data re-elaboration [42] is shown 
in figure 32 by blue squares. The global profile of 
the empirical <ln(A)> in figure 32 would resem- 
ble to that measured by Auger in a significant 
aspect: <ln(A)> decreases above the energy of 
10 17 eV up to (2 - 3) x fO 18 eV, where a mini- 
mum is attained. In the interval 2 x I0 18 -2 x I0 19 
eV, with some imagination, an increasing trend 
would faintly delineate, being in accord with the 
<ln(A)> profile of the Auger experiment above 
3 x 10 18 eV (see fig. 27). 

The partition of atmospheric showers in two 
groups according to the quoted re-analysis of the 
Agasa data [25] has been here refined with 6 par- 
ticles (or 6 ion groups) giving the outcome shown 
in figure 32 by pink square data points. In this re- 
finement the Theory of Constant Spectral Indices 
with the HE4 blend has been used. The resulting 
<ln(A)> profile (pink squares in fig. 32) is sys- 
tematically shifted downward by about one unit 



Figure 34. Measurements of the chemical com- 
position of the cosmic radiation by HiRes (pink 
dots) and Auger (black dots) experiments using 
the mean width of the longitudinal profile of the 
fluorescence light generated in atmospheric show- 
ers. The theoretical mean width in g/cm 2 versus 
energy evaluated by Heck [30] for protons and Fe 
nuclei with QGSjet and Sibyll algorithms is also 
shown as turquoise bands. 

of <ln(A)>. To appreciate numerically such a 
shift, just caused by the refinement with 6 ion 
groups instead of 2, consider that a Fe fraction 
of 0.60 implies a <ln(A)> of 2.41 (A = 56). The 
HEA blend at the energy of 10 19 eV for the 6 ions 
H, He, CNO, Ne-S, Ca and Fe with fractions, re- 
spectively, of .22, .21, .16, .16, 0.02 and 0.23 gives 
a <ln(A)> of 2.26. The LE blend at 10 19 eV with 
the ion fractions of .08, .4L .10, .09, 0.01 and 0.31 
gives a <ln{A)> of 2.43. U 

3 Unfortunately, an heavy ion fraction Fp e higher than 1 
results in figure 4 of ref. I42| , which signals a logical inco- 
herence whatsoever in the data elaboration materializing 
an unphysical situation. The incoherence is not easy to be 
eliminated since reverberates elsewhere 1461 . Probably, an 
additional correction of a few g / cm 2 in the Xf n f- X derived 
from the Mocca-Sibyll code would eliminate the Fe frac- 
tions higher than 1 in fig. 4 of the quoted reference I42| or 
in fig. 4 of ref. 44 . Note that this additional correction 
probably is unnecessary since the oversimplified data elab- 
oration with only two ion groups (instead of 6 or more) 
entails a bias in <ln(A)> which is clearly demonstrated 
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Figure 35. Evidence for the present uncertain- 
ties in cosmic-ray flux measurements in different 
experiments using giant terrestrial cascades in air 
above 10 17 eV. The spectrum (black curve) result- 
ing from the Theory of Constant Spectral Indices 
is normalized at 10 12 eV (black square) with a 
flux of 3.79 x 10 17 particles/(m 2 sr s eV 1 - 5 ) and 
it joins the Kaskade data in an excellent accord. 

A consistency verification of the chemical com- 
positions derived from the Agasa-Akeno data [43] 
and that extracted by X max measured by Fly's 
Eye has been performed in 1998 [44]. Both the 
radial muon density from the cascade core and 
the X^l x have been simulated with the same 
hadronic code, Mocca-Sibyll. Calculations have 
been adapted to Akeno (Al), Agasa (A100) and 
Fly' s Eye detectors including muon thresholds, 
trigger efficiencies and other instrumental con- 
straints, which are significant. Grouping all nu- 
clei in two fractions, light and heavy, so that F p + 
Fp e =l, it results a fraction Fp e greater than 1 for 
the Fly's Eye data around 10 17 eV. In order to re- 
move this incoherence a shift of - 30 g/cm 2 in the 
Xmax 1S required. The interesting aspect of this 
correction is that the theoretical profile X^ x re- 
sulting from Mocca-Sibyll code differs from that 
used in this paper (see figure 12) by about 25 



in figure 32 or in figure 33 with the Volcano Ranch data. 



g/cm at 10 eV, and similarly at higher en- 
ergies. Such a difference would alleviate, as a 
gratuity, the gap between the <ln(A)> predicted 
by the Theory of Constant Spectral Indices with 
r=2 and the Auger data (see fig. 27). 



An elaboration and revision of atmospheric 
shower data of Volcano Ranch experiment oper- 
ated in the period 1959-1964 has been performed 
in 1975 using a bi-modal analysis [35] . Figure 
33 shows the theoretical profile (red curve) along 
with the data point of Volcano Ranch (blue rect- 
angle) spanning the interval 5 x 10 17 -10 19 eV. The 
resulting <ln(A)> has a mean value of 3, which 
indicates that the cosmic radiation at 10 18 eV is 
still dominated by intermediate and heavy nu- 
clei. If the same procedure of ion ponderation 
applied to the Agasa data is also applied to Vol- 
cano Ranch data with the HEA ion blend, the 
mean value of <ln(A)> decreases from 3 to 2.4 
(pink cross in the rectangle) approaching the the- 
oretical value of 2.5 as displayed in figure 33. 



10. Reported measurements of the mean 
atmospheric depth of the cosmic radi- 
ation by fluorescence light in HiRes 



Though the longitudinal profile of giant atmo- 
spheric showers reconstructed by the fluorescence 
light is rather insensitive to the hadronic mod- 
els employed in the data analysis, it would seem 
that some instrumental unknowns or uncontrolled 
measurement procedures still plagues the deter- 
mination of X max by fluorescence light. 

Preliminarily, notice that experiments taking 
advantage of the fluorescence light in the mea- 
surements of X max like HiRes do not agree one 
another. The X max observed by the HiRes exper- 
iment [38] is larger than that measured by Auger 
[TB] which uses the fluorescence light as well. As 
an example, the Fly' s Eye data on <ln(^4)> [3"T] 
in the energy decade 3 x 10 17 -3 x 10 18 eV shown 
in figure 28 differ from those of HiRes Prototype 



(figure 29) by one unit of <ln(A)> which is quite 
large by any standard, theoretical, empirical or 
instrumental. 

The measurements of X max by florescence light 
demands the interpolation of the longitudinal 
profile by an appropriate function fop which 
has a characteristic width denoted here u{X max ). 
The average value of the distribution of a(X max ) 
at a given energy is a measurement of the chem- 
ical composition. The observed u(X max ) may be 
cross-check method in the measurements of the 
chemical composition with the same data sam- 
ples. 

Presently (2009), the available measurements 
of a(X max ) of the HiRes Collaboration [47] dis- 
agree with the analogous data of the Auger ex- 
periment 48 as shown in figure 34. Note that 
the theoretical <j{X max ) profiles shown in figure 
34 derive from the same sample of the simulated 
cascades used to calculate the X max profiles of fig- 
ure 12. The (j{X max ) profile measured by HiRes 
is rather flat around the value of 50 g/cm 2 in the 
interval 10 18 - 5 x 10 19 eV while Auger data (full 
black dots) exhibit a clear decreasing profile from 
a maximum of 58 g/cm 2 down to a minimum of 
22 g/cm 2 . Imagine a linear scale of energy be- 
tween 10 18 - 5 x 10 19 eV in figure 34 and the dis- 
agreement between the two experiments dilates, 
since the initial agreement below 4 x 10 18 eV is 
confined in less than 5 per cent of the explored 
energy band. 

At the energy of 10 17 eV the profile of <ln(A)> 
measured by Kascade joins that of the HiRes ex- 
periment [35] (fig- 29). The average value of 
the <ln(A)> of the Kascade experiment is about 
3.1 (average value of 3 methods of measurements) 
while that of HiRes is close to 2.5. It is important 
to emphasize the disagreement, or the potential 
disagreement, between these two experiments on 
<ln(A)> around 10 17 eV. The knees of the light 
ions have been observed by a number of exper- 
iments above 10 15 eV [4=9)50] . This fact neces- 
sarily implies an increase of <ln(A)> above 10 15 



eV, since the disappearance of light ions automat- 
ically enhances the heavy ion fraction. Knowing 
that at 10 15 eV the value of <ln(A)> is certainly 
comprised between 1.74 and 1.78 as inferred by 
extrapolating balloon and satellite data [52], any 
estimate of <ln(A)> at 10 17 eV leads to the value 
of 3.1 and not a lower value as HiRes Prototype 
data (2004) [36] around 10 17 eV in figure 29 would 
suggest. 

As apparent in figure 29 the mean values of 
<ln(A)> of the HiRes experiment after the ac- 
ceptance correction [38] are globally displaced 
downward by some 0.5 units. The shift is surpris- 
ingly large compared to the reported error bars 
before and after acceptance corrections [38] . The 
large excursions in <ln(A)> (blue dots) also sig- 
nal that the global detector performance suffered 
from residual uncontrolled unknowns. 

The correction of the detector acceptance [33] , 
a routine effect applied after 3 years from the orig- 
inal measurement [37] does not preclude future 
data re-analysis. 

In synthesis, the HiRes data pertaining the 
chemical composition of the cosmic radiation dis- 
agree 0: (1) with the Auger data on X max ; (2) 
with the Auger data on cr{X max ) (see figure 34 
and the related text); (3) with those of all other 
experiments on <ln(A)> above 10 17 eV (Volcano 
Ranch, Haverah Park, Akeno, Agasa, Fly' s Eye 
and Yakutsk); (4) with the Kascade data at 10 17 
eV. Theoretically, a bias exists in the conversion 
of X max into <ln(A)> made by the HiRes Col- 
laboration because, out of all current versions of 



4 The resonance of this conclusion tunes up with the crit- 
icism expressed by A. A. Watson I51| on Fly's Eye and 
HiRes experiments: <C-... By contrast, with florescence 
detectors the aperture continues to grow with energy and 
remains considerable uncertainty about the HiRes aper- 
ture. Extensive Monte Carlo calculations must be made 
to establish it and these make assumptions about the slope 
of the spectrum and about the primary mass, with pro- 
tons being assumed on the basis of the claims from the 
Fly' Eye and HiRes experiments ignoring often contradic- 
tory evidence from other experiments without discussion 
>. 



hadronic codes adopted to simulate nuclear in- 
teractions in air (e.g. QGSjet, Sibyll, etc.), those 
adopted by HiRes tend to generate protons and 
depress heavy ions (see fig. 8 of ref. 10 and fig. 4 
ref. [39]). 



e.g. r=I ed /I d . 

From the ensemble of the profiles of <ln(A)> 
extracted from the measurements of X max of 
the Auger, Yakutsk, Fly' s Eye, Agasa, Akeno, 
Haverah Park and Volcano Range experiments it 
emerges a global accord with the theory. 
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Figure 36. Theoretical difference of the atmo- 
spheric depth for protons and Fe nuclei according 
to a number of evaluations. The model QGSjet-03 
(thick pink curve) would reduce the discrep- 
ancy between Auger data on <ln(A)> and the 
corresponding profile derived from the Theory of 
Constant Indices. 

11. Conclusions. 



The comparison of the chemical composition of 
the cosmic radiation expressed by <ln{A)> de- 
rived from the Theory of Constant Indices 
with the observed <ln(A)> extracted from the 
Auger data on X max dictates the necessity of in- 
troducing an extradisc component of the cosmic 
radiation above 10 17 eV. In Sections 1 and 2 it 
is explained why this component is better termed 
extradisc component, I ed . The variable r is a free 
parameter of the theory defined as the extradisc- 
to-disc flux ratio in the solar cavity at 10 19 eV, 



In the energy interval 10 15 -10 17 eV there is an 
excellent agreement with the observations of Eas- 
top, Kascade, Tunka (see fig. 3) and other exper- 
iments. Below 10 17 eV the extradisc component 
is a negligible fraction of the disc component Id, 
and as a consequence, the theoretical and empir- 
ical analysis reported in the companion paper [8] 
remains valid. 

Because of the unknowns entailed in the con- 
version of X max into <ln(A)> above 10 17 eV, the 
accord between the theoretical profile and data in 
figures 27, 28, 30, 31, 32 and 33 is more than satis- 
factory. Values of r= I e d/Id higher than 2 would 
alleviate the gap between the empirical profile of 
<ln(A)> and the theory, as shown in figure 27. 
Presently, values of r higher than 2 are neither 
convincing nor necessary, due to large systematic 
errors in the experimental data. The last two fig- 
ures 35 and 36 intend to illustrate this position. 

Figure 35 depicts in a flash the present experi- 
mental uncertainties in flux measurements quot- 
ing arbitrarily data from Kascade, Yakutsk and 
Auger experiments. Fluxes from other experi- 
ments do not agree as well. Such uncertainties 
impede to identify a more precise and reliable 
value of the ratio r in figure 26 on an empiri- 
cal basis, being r=2 the present adequate global 
estimate. 

Let be X™ l ax (H) the theoretical atmospheric 
depth for protons, X™ l ax (Fe) that for Fe nuclei 
and D max the difference between the two depths 
i.e. D max = X^ ax (H)- X^ x (Fe) at a given 
energy E. Error sources generating differences 
in <ln(A)> in theories and in experiments might 



arise both from X^ x (H) and D max as well. Ide- 
ally, different haclronic models may result in a 
combination of X^ x (H) and D max having the 
same <ln(A)>. Figure 36 reports the D max ver- 
sus energy in some current hadronic models and, 
at the arbitrary energy of 10 175 eV, the X^ x 
(H). Figures 36 and 12 show that the X%* x (H) 
and D max profiles used in this paper for the con- 
version of X max into <ln(A)> are quite similar 
to those adopted by others (for example, Auger). 

Notice that the hadronic code called QGSjet 
model-03 [52], having a larger D max of 120-130 
g/cm 2 in the relevant range 10 17 -5 x 10 19 eV and 
a deeper X™ l ax (H) of 739 g/cm 2 at 10 175 eV (see 
figure 36), it would shift upward the simulated 
Xmax profiles in figure 12 by 9 g/cm 2 at 10 17 ' 5 
eV. The same correction applied to the Auger 
data in figure 3 (or figure 14) would shift upward 
by 0.57 units of <ln(A)> (9 g/cm 2 for r=0 or 
11 g/cm 2 for r=2) at 4.12 x 10 17 eV alleviating 
the gap between Auger data and theory. As far 
as uncertainties in hadronic models above 10 17 
eV used in the analysis of giant atmospheric cas- 
cades remain large, the determination of a value 
of r higher than 2 based on the <ln(A)> data 
would appear premature. This same conclusion 
is suggested by the flux measurements in figure 
35. 

Let us finally mention the systematic disagree- 
ment between the <ln(A)> extracted from the 
Xmax measured by the HiRes experiment and 
theory, amounting to about one unit in the en- 
tire interval 10 17 -5 x 10 19 eV (see fig. 29). The 
examination of the HiRes data on <ln(A)> made 
in Section 9 suggests a number of inconsistencies 
with the data of numerous experiments and with 
different versions of the HiRes experiment itself 
(HiRes Prototype, HiRes Stereo, data revision 
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